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Abstract 

The intrinsic decoherence from vibrational coupling of the ions in the Cirac- 

Zoller quantum computer [Phys. Rev. Lett. 74, 4091 (1995)] is considered. 

Starting from a state in which the vibrational modes are at a temperature 

T, and each ion is in a superposition of an excited and a ground state, an 

adiabatic approximation is used to find the inclusive probability P(t) for the 

ions to evolve as they would without the vibrations, and for the vibrational 

modes to evolve into any final state. An analytic form is found for P(t) at 

T = 0, and the decoherence time is found for all T. The decoherence is found 

to be quite small, even for 1000 ions. 

PACS numbers: 89.80.+h, 32.80.Pj, 03.65.Db, 36.40.+d 
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Quantum computers (QC) are (as yet hypothetical) devices with states that are quan- 
tal in nature, and which perform calculations by unitary transformations on these states 
The linearity of the superposition principle leads to an inbuilt massive parallelism: a 
computer with N two-state elements can operate on 2^ states simultaneously. This paral- 
lelism underlies Shor's recent algorithm || for factorizing a composite number of order 2 L 
in ~ L 3 steps on a QC. The best known classical algorithm takes ~ exp[c(L)L 1//3 ] steps with 
c(L) ~ (InL) 2 / 3 . The potential for other quantum algorithms is clearly exciting, as is the 
emergence of a new paradigm for computation itself. 

It is obvious that maintaining perfect phase coherence among all the states of a QC is a 
daunting task, not to mention getting these states to evolve in the desired fashion in the first 
place ||. An imaginative proposal for a QC by Cirac and Zoller || (CZ) seems promising in 
addressing these problems [0|§. It utilizes a string of N identical ions in a linear Paul trap 
@, with each ion separately addressable by a laser. Two internal states of each ion, |e) and 
| g), are used for the QC, along with the center-of-mass (CM) axial vibrational mode of the 
entire array. A program is implemented as a specified sequence of (vr/2, ir, etc.) pulses that 
drive |e) «-> \g) transitions on any given ion, along with pulses detuned by the CM frequency 
that enable coupled transitions between any pair of ions. 

Two types of decoherence should be distinguished in the CZ (indeed, any) QC. The first 
is technical, due, e.g., to imperfect phase locking, mistuning of lasers, errors in timing and 
duration of pulses, and overlooked perturbations in the Hamiltonian. The second kind is 
intrinsic, and arises from coupling of the computationally useful to the undesirable bath 
degrees of freedom. Although the technical problems alone render the pursuit of a QC a 
fool's quest in many people's eyes, intrinsic decoherence sets basic limits on the capabilties 
of a QC. It is with this motivation that we study intrinsic decoherence in the CZ QC fllPf . 

We take as our bath the vibrations of the ions, which we treat as undamped harmonic 
oscillators. Damping can be included if necessary ]TI[. Radiative decoherence is accounted 
for very simply by demanding that any computation take less time than r sp /N, where r sp 
is the spontaneous |e) — > \g) decay time for one ion. It clearly pays to have as large a r sp 
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as possible, by working with El forbidden transitions ||, or with hyperfine sublevels of the 
ground ionic multiplet 0. The total Hamiltonian minus the driving lasers can be generally 
written as (setting % — 1) 



1 p 2 1 

Here, the cr's are equivalent Pauli spin operators in the {|e), \g)} space, a± = (a x , a y , 0), 
and p^ are the vibrational normal mode coordinates and momenta, and m is the mass of each 
ion. We shall refer to the three terms in Eq. ([!]) as Hi, H nm , and H' respectively. The Cj M 
are calculable functions (see below) of the ionic transition matrix elements and equilibrium 
ion positions, which we assume are such that there is no \e) <-> \g) transition term in H in 
equilibrium. It is key to successful operation of the CZ QC that the vibrations be cooled to 
nearly zero temperature, and that the frequencies and the couplings c ifJL be small. The 
approximations of this paper require that u ^> u; M , (X^Q/^), and J2 u{ c in%) 2 1 luj o oj^ 
which as we shall see, can be satisfied comfortably. 

Let us now study the effects of the bath on the simplest computation of all, i.e., just 
waiting. We include the CM vibrational mode in the bath for simplicity in this note, as 
this is not expected to change the result qualitatively. Suppose that initially, the bath is 
described by a density matrix p(Q, Q') (Q denotes all the collectively), and the ions are in 
some state |in). The system is not driven by any lasers, and simply sits for a time t. What 
is the probability P(t) of finding the ions in the final internal state |fin) = exp(— iHit)\m) 
that one would get in the absence of the bath, and the bath in any state whatsoever? For 
the state |in), we take 

|m> = Il2- 1/2 (l+>* + !->,), (2) 

i 

where a iz \±)i = ±|±)«. This state is illustrative of the complex superpositions of computa- 
tional basis states that give QC's their parallelism. Since Eq. (|I|) describes a finite, closed 
system, P(t) -/-> as t — ► oo, but we expect that P(t) will drop close to zero at some time 
Ta, after which it will fluctuate with small amplitude [|l^]. The time limits the longest 



computation that can be done with the CZ QC (if t c i < T ap /N). The coherence time is 
expected to decrease when transitions are driven by the lasers, and can also be estimated 

To evaluate P(t), we write the reduced bath density matrix propagator as a double path 
integral 

rQ'f rQf 

J(Q f ,Q' f ;Q,Q')= / 7 [dg][dQl e i (*WW]-*WW]U[g(t)]A*[g / (t)], (3) 

JQ> JQ 

where Sq is the action for the bath alone, and 



Texp (-ij [Hi + H'(Q(t'))} dt' 



in)- (4) 



A[Q(t)} = ^fin 

In terms of J, P(t) is given by 

P(t) = JJJdQjdQ dQ' J(Q f , Q f ; Q, Q') p(Q, Q'). (5) 

It now does no good to integrate out the oscillators. Instead, we exploit the fact that 
u ^> uj^ — typically u ~ 10 15 s _1 , and ~ 10 7 s _1 — to integrate out the spins. It is 
easily seen that factorizes into Yii Ai[Q(t)]- The zth spin experiences a field u>ix_ = 

2J2fj, c i^<l^(t) m the x-y plane, with magnitude Ui± <^ ujq and \dlnu>i±/dt\ u>o, as we may 
safely assume that Eq. (§) is dominated by paths Q(t) varying slowly on the Uq 1 time scale. 
This permits us to evaluate v4j[Q(t)] using an adiabatic approximation. We may further 
take the instantaneous precession axis of the spin as z for all t with negligible error (of 
order to^/tOo). It is far more important to approximate the phase well. The instantaneous 
energies of the states |±)j are given by ±(u + ufj_/2u ) to relative order Ui±/u . We thus 
obtain Ai = cos$j(t) with 

fdtulJfi/Auo. (6) 
Jo 

Equation (0) holds for t < O (8uJq / (uf ± )) , where () denotes an average value, and the dot 
denotes d/dt. (The results obtained below imply that this time scale exceeds r^.) We have 
also omitted a negligible Berry phase term |TT| . 



It now pays to rearrange the expression for P(t). First, we write cos $j as a sum of e ±l *% 
and substitute the resulting expression for A[Q(t)] in Eq. (|3|). This yields, 

J(Q f ,Q> f -Q,Q>) = -L J2 K {s} {Q f ,Q)K\ s , } {Q' f ,Q'), (7) 

with 

K{ s }(Qf,Q) = f Q, [dQ] exp liS [Q(t)\ +i £ ^ f q,(t)q u (t) dt\ . (8) 

In Eqs. (0) and (|8|), {s} = (si, s 2 , . . . , Sjv), {s'} is similarly given, each Sj = ±1 is an Ising-like 
variable, and u* = c ilM -c iu /u . Next we define the following combination of propagators, 

R {s>sl} (Q,Q') = J dQ f K {s} (Q f ,Q)Kl sl} (Q f ,Q'), (9) 

in terms of which P(t) can be written as [see Eqs. @ and (|7])] 

p W = if E II dQdQ'R {S)Sl} {Q,Q')p{Q,Q'). (10) 
2 {«},{«'} 

It is apparent that -ft'js} is the propagator for a set of coupled harmonic oscillators, 
described by a Hamiltonian that depends on the Ising congiguration {s}: 

hu{s}) = E ll + 1 E <^ 2 (ii) 

(fiV = ^-2E s ^;,- (i 2 ) 

i 

We can thus write R{ s , s '} alternatively as 

R{s,s'}(Q,Q') = (Q'le^^e-^^IQ). (13) 

Since Eq. (|l"3"|) only involves harmonic oscillators, we can evaluate it exactly by reverting 
to path integrals. The exact answer involves trigonometric functions and determinants of 
the matrices fi({s}) and fi({s'}) and is of limited use because of the remaining sum on the 
s's. To make further progress, we employ an approximation in the same spirit as that used 
to obtain the spin transition amplitude Aj(Q[t]). Namely, accuracy in the normal mode 
frequencies is much more important than in the normal modes themselves. Errors in the 



former lead to errors in and R{ s , s '} that grow with time, while errors in the latter do not. 
We therefore treat the second term in Eq. fll2| ) as a perturbation, and use the unperturbed 
normal modes, but correct the frequencies to first order: oj'^ = — J2i s i ^LtAv [Given 
the stated assumptions about the relative sizes of the three terms in Eq. ([1]), the frequency 
shift can indeed be seen to be small.] With this approximation, the kernel R{ s , s >} factorizes 
into a product Rp of kernels for each mode (we suppress the Ising variables where no 
confusion is possible), with 

= (^^'"V^^lg,), (14) 

and Hp = p 2 j2m + muj'^q^/2. Equation (|T4]) has a simple physical interpretation. Starting 
from an initial state, the system evolves forward in time for a duration t as a harmonic 
oscillator of frequency iO\ , say. It then evolves backward in time for duration t as a harmonic 
oscillator of slightly different frequency u>2. For our problem, this difference propagator, R^, 



can be further simplified because the frequencies U\ and u 2 are almost identical []T3|]. If we 
think about the corresponding classical problem in phase space, the forward and backward 
evolutions take place on ellipses of nearly equal eccentricity. To good approximation, we may 
regard the ellipses as coincident. With suitably scaled p and q axes, this common ellipse is 
a circle, on which the particle sweeps out angles uj\t and — u 2 t in the forward and backward 
motion. The net evolution is that of a single harmonic oscillator of frequency 8^ — u 2 — &>i, 



and mass for a time t fL4j . In other words, 

RM^%) « ( g ;| e ^+™^)V/2m"M| gM ). (15) 



It is now easy to carry out the coordinate integrals in Eq. ( fTCf ) for the special case 
where p is a thermal equilibrium density matrix oc e~@ Hnin with j3 = 1/kT. Since p and 
R both factorize by normal mode, i.e., p = EL/^, and R = n^R^, the summand in 
Eq. (p!0|) also factorizes into H^A^t), where A M = // R^p^dq^dq^. By using standard 
coordinate representations of the harmonic oscillator density matrix and propagator, is 
easily evaluated, and the result can be written as (restoring %) 



P(t) = 2-™ £ l[A,(t), (16) 

{s,s'} A* 

A m = sinh/^/2 

M ; sinh(/?^-?V)/2' 1 J 

Note that 5^ = Y,i( s i ~ s 'd ul u,J u n depends on the Ising configuration, and u* is given below 
Eq. (|). 

Equations (|T6|) and (0) formally answer the question we set out to investigate, but the 
sum on the s's is nontrivial. Some general properties of the result are worth noting, however. 
Thus, P(t) is real, and since |A M (i)| < 1, P(t) < 1. [In fact P(t) = 1 only if i/m^ 
is an integer for all /i simultaneously] The expressions simplify greatly at T = 0. Then, 
A^(£) = e ^*/ 2 , and 

P(t)=ncos 2 (C 2 t/2^ ), (T = 0), (18) 

i 

(i = E C = E <v- <v/mnay (19) 

Note that £j and have dimensions of frequency. We can also obtain the decoherence time 
Td for all T by examining the initial drop of P(t) from unity. Writing 1 — P(t) ~ (i/r^) 2 , we 
obtain 
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4u; 2 ^ 



\ 2 

E C> coth xA + E Ccosech 2 ^^ 



(20) 



i L v 7 fi 

where = (3fiuj^/2. Note that falls as T rises, as it should. At T = 0, tJ 2 = X^CV^o- 
We still need the couplings c ifl . These depend on the nature of states |e) and \g), so we will 
find them only for a particularly favorable situation obtained by using Ba + ions, and states 
with AM = ±1 in the 6s 2 S\/2 and 5d 2 D§/ 2 (or 2 -D 3 / 2 ) multiplets for \g) and |e) respectively. 
The 2 D 5 / 2 — > 2 5'i/2 decay is an E2 process, with r sp ~ 35 s |15|1 , co> /27r = 1.7 x 10 14 Hz. The 
interaction hamiltonian is given by 

h ' = qY! E % 4 <^L, (2i) 

i,j a=x,y 

where g is the ionic charge, z iy Q l a/3 , and u l a are the equilibrium position, quadrupole moment 
tensor, and displacement from equilibrium of the zth ion; Zij = Zi — Zj, u % i = u\ — u J a . Note 



that only transverse vibrations appear in Eq. ( |21| ) because states |e) and \g) are connected 
by | AM | = 1; the same restriction ensures that the equilibrium quadrupole fields do not 
drive any |e) <-» \g) transitions. Writing the mode index \x = (r, a), with r = 1, 2, . . . , N, we 
can write = E r F^q ra , where the F r * are normal mode eigenvectors |J16[. We scale these to 
obey E r K F r = 8*, Ei^W = 5 rs . With Fj? = F* - F>, and Q l a/3 = (e|Q l J<?>, we obtain 

| c i,ra| = ~QY2j F r J z ij \Qaz\- 

It remains to substitute the above expression for Cj >m into Eqs. (|19D and ( f20|) to obtain 
Td- Here we will only outline a simplified calculation at T = 0. The steps are as follows. 
First, in Eq. flUf) we replace ^ by a^jy, the frequency of the zigzag transverse mode, for all 
fi. Since uj t .N is the smallest normal mode frequency, this replacement yields a lower bound 
for Td- It also obviates finding the F*, as the sum over r can be done by orthonormality 
of the F r \ Second, for the states in question |<3^J 2 = 18^Ao r sp> with k = uo/c, and 
a — x, y. Third, we note that if u z is the longitudinal CM vibrational frequency, d = 
(g 2 /™^) 1 ^ 3 is a natural trap length scale. We define the dimensionless sums S n (i) = 
Y,'j z ij n dQ. [The expression for t^ 2 contains a factor E« (^IW + ^(z)) 2 .] We estimate S n (i) 
using a continuum approximation for the ion array |I7|. In this approximation, the local 
interionic spacing is s(z) = s(0)(l — 2 2 /L 2 ) -1 , where s(0) = 4L/3N, and 2L is the total 
length of the array, with L 3 « 3iVln(0.8iV). (All lengths are in units of d .) This yields 
S n (i) ~ 2((n)/s n (zi). The sum over i can now be estimated by an integral. Combining 
these results, we obtain a bound for r d entirely in terms of trap and ion parameters, 
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1 1 iV 35/6 
— < n qfi z (22) 

r d ■ T sp [\n(0.8N)}^cu uj t , N (W 1 ' 

The last step is to estimate oot,N- By considering the transverse force on the central ion, 
we get uj 2 n ~ cu 2 — c(iV)u; 2 , where cjj is the CM transverse vibrational frequency, and 
c(N) = 9C(3)jV 2 /16 ln(0.8iV). In fact, requiring u t ,N > gives the critical u t {oJ t>CT ) needed 
to avoid the zigzag instability [|T7]]. The numerical value of r d implied by Eq. ( p2]) is a 
very sensitive function of u z (~ oo~ ~ 16 / 3 ). By choosing u z /2tt in the 10-100 kHz range, and 
uo t) N — ^t,cr/2, the ratio Td/r sp can be seen to lie in the 10 4 -10 8 range for N = 1000. [The 
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assumptions behind Eqs. (|j), fll4|), and (|15|), can all be seen to hold.] 

We thus see that contrary to what might have been expected, vibrations of the ions are 
not a significant source of decoherence in the CZ QC for iV < 1000. Larger iV values pose 
serious technical challenges in trap design and in keeping s(0) large enough to be optically 
resolvable. The conclusions of the present paper, however, can only be encouraging for the 
prospect of quantum computation. 

I am grateful to H. J. Kimble and D. Wineland for useful comments and correspondence. 
This work is supported by NSF grant no. DMR-9306947. 
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